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ABSTRACT 

GRB afterglows offer a probe of the intergalactic medium out to high redshift which 
complements observations along more abundant quasar lines-of-sight. Although both 
quasars and GRB afterglows should provide a-priori random sight-lines through the 
intervening IGM, it has been observed that strong Mg-II absorbers are twice as likely to 
be found along sight-lines toward GRBs. Several proposals to reconcile this discrepancy 
have been put forward, but none has been found sufficient to explain the magnitude 
of the effect. In this paper we estimate the effect of gravitational lensing by galaxies 
and their surrounding mass distributions on the statistics of Mg-II absorption. We 
find that the multi-band magnification bias could be very strong in the spectroscopic 
GRB afterglow population and that gravitational lensing can explain the discrepancy 
in density of absorbers, for plausibly steep luminosity functions. The model makes 
the prediction that approximately 20%-60% of the spectroscopic afterglow sample (i.e. 
^ 5— 15 of 26 sources) would have been multiply imaged, and hence result in repeating 
bursts. We show that despite this large lensing fraction it is likely that none would 
yet have been identified by chance owing to the finite sky coverage of GRB searches. 
We predict that continued optical monitoring of the bright GRB afterglow locations 
in the months and years following the initial decay would lead to identification of 
lensed GRB afterglows. A confirmation of the lensing hypothesis would allow us to 
constrain the GRB luminosity function down to otherwise inaccessibly faint levels, 
with potential consequences for GRB models. 

Key words: cosmology: theory - gamma rays: bursts -quasars: absorption lines - 
gravitational lensing 



1 INTRODUCTION 



for cross-correlation yields a statistical accuracy of a factor 
of two in host halo mass. 



At moderate redshifts (0.3 < z < 2.2) the dispersal of metal 
enriched gas through the IGM via galactic winds is most 
readily traced via Mg-II absorption systems in optical spec- 
tra of bright background sources. Although the detailed ori- 
gin of Mg-II absorbers remains uncertain it is thought that 
they are associated with galaxies and galactic outflows. For 
example, Mg-II absorbers have been shown to be associated 
with neutral hydrogen absorbers over a range of column 



densities, including damped Lyman-a absorbers (e.g. Rao 
et al.|2006 l. Moreover the host halo masses associated with 
Mgll absorbers have been estimated at 2 ~ 0.5 via cross- 
correlation with luminous red galaxies in the Sloan Digital 
Sky Survey Data Release 3 ( Bouche et al.|2006 l, yielding a 
host mass of M ~ lO^'^M© (i.e. massive galaxies). Indeed 
the large number of Mgll absorbers and galaxies available 



Obtaining an unbiased census of the density and distri- 
bution of Mg-II in the IGM requires that the back-ground 
sources be uncorrelated with the foreground absorbers under 
study. Most current knowledge regarding the census of Mg- 
II absorption systems comes from spectra of quasars (e.g. 
Prochter et al. 20061. More recently, observations of GRB 



afterglows have begun to offer a new probe of the intergalac- 
tic medium out to high redshift which complements the more 
abundant quasar lines-of-sight. Indeed, since they are asso- 
ciated with star formation rather than supermassive black 
holes (which, given their long assembly times, become ex- 
tremely rare at high redshift), GRBs could potentially be 
seen out to much higher redshift than quasars; the current 
record holder is at 2: ~ 8.1 ( Salvaterra||2009 1 . Interestingly, 
although both quasars and GRB afterglows should provide 
a-priori random sight-lines through the intervening IGM it 
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has been observed that strong Mg-II absorbers are several 
times as hkely to be found along sight-lines to GRBs as 



of high resolution spectroscopy from Prochter et al. ( 2006 1 



along quasar sightlines (Prochter et al. 2006 Tejos et al 



2009 Vergani et al. 20091, indicating that one or both of 



these samples are biased relative to the underlying popu- 
lation of absorbers. Several proposals to reconcile this dis- 
crepancy have been put forward, with a detailed discussion 
of possible biases presented by Porciani et al. (20071. For 



example, the incidence of Mg-II systems in quasars could be 
lowered because of dust obscuration associated with the ab- 
sorbing systems; this turns out to be too small to explain the 
discrepancy. Alternatively, it has been argued that different 
sizes of the source could lead to different absorber incidence 



between GRB afterglows and quasars (Frank et al. 20071 



However, from the similarity of the equivalent width distri- 



butions in GRBs and quasars, Porciani et al. (20071 show 



that the absorbers must be larger than the beamsize, and 
hence it is not possible to explain the difference in this way. 

A third potential bias is provided through gravitational 
lensing of GRBs by foreground galaxies associated with the 
Mg-II absorber. There are several lines of circumstantial evi- 
dence for this. Firstly, imaging studies of the host galaxies of 
GRBs with early time afterglow spectra ( Chen et al.||2009 l 
show that additional galaxies are found within 2 of all GRB 
host fields where the line-of-sight contains a Mg-II absorber, 
while no additional galaxies are found within 2" of GRB 
lines-of-sight that do not have any Mg-II absorbers. More- 
over, GRB afterglows that have more than one absorber are 



found to be a factor of 1.7 brighter than the others ( Porciani 
et al.|2007 1 . Under the assumption that gravitational lensing 



magnification of GRBs is responsible for the excess number 
of absorbers in GRB spectra relative to quasars, |Tejos et al.| 
(20091 estimated the fraction of GRB afterglows that need 



to have been included in the sample owing to magnifica- 
tion. They find a value of 60% which indicates a large grav- 
itational lensing induced magnification bias. On the other 
hand, no strongly lensed GRBs have been discovered]^ At 
first sight this suggests that gravitational lensing could not 
provide the required explanation for the excess absorbers in 
GRB afterglow spectra. However, as noted by [Porciani et al^ 
(20071 the selection of GRB afterglows in two (apparently) 
independent bands (namely gamma-rays and the optical) 
leads to the possibility of an increased multi-band magnifi- 



cation bias (Wyithe et al. 20031. One goal of this paper is 



to determine under which circumstances a sufficiently large 
magnification bias can be obtained. 

The most recent collection of GRB afterglow spectra 
from which strong (> lA) Mg-II absorption can be stud- 



ied is summarised in the work of Vergani et al. ( 2009 1 who 



compiled a list of 26 spectra (with a combined redshift path- 
length of Az = 31.55) containing 22 strong Mg-II absorbers 
among 15 of the 26 lines-of-sight. From that paper we take 
the following values for observables describing the absorber 
population. The fraction of lines-of-sight that contain one or 
more absorbers is ~ 0.6±0.15. By contrast, the results 



^ The expected rate of strongly lensed GRBs was discussed by 
[Porciani &: Madau| | [2001[ | , who found that although a few strongly 
lensed GRBs should be present in a 3 year survey with Swift, 
the partial sky coverage means that it is unlikely that lensed pairs 
would have been identified as recurrent bursts. 



yielded 22 absorbers along 91 quasar lines-of-sight, implying 
that Jivig.q = 0.25. Assuming the same redshift pathlength 
distribution as the GRB sample, these values imply that the 
relative incidence of the number of strong systems in GRB 
afterglows relative to quasars is -RoRB.q ~ 2.1 ± 0.6. Where 
required we take the redshifts of GRBs and associated ab- 
sorbers, as well as the probed redshift path-lengths from this 
paper. Our goal is to make a quantitative comparison of a 
simple parameterised model for the expected number of ab- 
sorbers (which includes the bias introduced by gravitational 
lensing), with these observations. 

Our paper is presented as follows. Section [2] describes 
the basics of the lensing model with which we interpret the 
absorber statistics. In §[3]and §|4]we discuss the statistics of 
Mg-II absorbers in GRB afterglow spectra, and their com- 
parison with statistics from quasar lines-of-sight within the 
context of our model parameters. We then discuss the pre- 
dicted distribution of absorber redshifts (§ |5| and observed 
source mag nifications (§ |6|. In § [t] we discuss the rate of 
multiple imaging among the spectroscopic GRB afterglow 
sample before presenting our summary in § [8] In our nu- 
merical examples, we adopt the standard set of cosmological 
parameters ( Komatsu et al.|2009 1, with values of fib = 0.04, 
Q.ra = 0.24 and Q.q = 0.76 for the matter, baryon, and dark 
energy fractional density respectively, and h = 0.73, for the 
dimensionless Hubble constant. 



2 LENSING MODEL 

The basis of our model is that we assume all lines-of-sight 
within a distance 



R = Ro 



(1 + ^)" 



(1) 



V200km/s^ 

of the nearest galaxy (e.g. Tinker &: Chen||2008 1 contain a 
Mg-II absorption system with an equivalent width greater 
than some specified minimum (taken as > lA for the Mg-II 
absorption statistics modelled in this paper), while lines-of- 
sight with impact parameters larger than R contain no ab- 
sorbers. This scaling is chosen to represent a constant mul- 
tiple of the virial radius. Here Ro is the proper radius within 
which a strong absorber will be found around a 200km s~^ 
halo at z = 0. We then evaluate the probability of a partic- 
ular line-of-sight containing a Mg-II absorber, including the 
effects of magnification bias owing to gravitational lens in- 
duced magnification of the observed source flux. For simplic- 
ity we model galaxies as singular isothermal spheres (SIS), 
which is an adequate approximation for the purposes of 
this paper ( jTreu fc Koopmans||2004| . In the case of close 
alignment, the presence of an intervening galaxy can pro- 
duce strongly lensed multiple images of background sources, 
which are magnified relative to their intrinsic flux (with a 
magnification fi > 2 for the bright image in the case of an 
SIS). Images that are not so closely aligned, and therefore 
only singly imaged can also be magnified (in the case of an 
SIS by a value of 1 < /i < 2). In this paper we distinguish 
between the numbers of absorbers in the spectra of multiply 
imaged and singly imaged sources. 

We refer to the a-priori probability for a source to be 
multiply imaged as the multiple image optical depth Tm- 
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In addition to this quantity we calculate the optical depth 
TMg for Mg-II absorption (i.e. the probability that a line- 
of-sight goes within i? of a galaxy), as well as the optical 
depth for which there are no absorbers TnoMg = 1 — TMg. 
The optical depth must be computed over the redshift path- 
lengths where Mg-II absorption has been searched. In the 
case of GRB afterglow spectra this is object dependent since 
different GRBs were observed with a range of telescopes and 
spectral ranges. Accounting for this, the multiple image and 
Mg-II absorption optical depths are 



(2) 



and 



TMg 



A''grb 



Ngrb 

E 



dz I da^a){l + zf'^nR'ia) 



(3) 

respectively, where 6'er is the Einstein radius and Dd is the 
angular diameter distance to the lens. Here we assume that 
the SDSS velocity function ( Sheth et al.|2003 \ describes the 
galaxy population (and assume a constant co-moving den- 
sity at fixed a) 



(4) 



where = 2 x 10"^Mpc"^, a = 6.5, P = 14.75/q and 
0-* = 16ir(Q/^)/r((a + l)//3). To calculate we use the 
expression for the angular Einstein radius for an SIS 



6Ier(o-) = 1.4' 



1)7 V200km/s 



(5) 



where Ds and Dds are the angular diameter distances to the 
source, and between the lens and source respectively. 

The GRB afterglow sample is selected from the wider 
GRB population using two flux limits. Firstly there is a se- 
lection in gamma-rays, and secondly a selection for optical 
afterglows that are bright enough to facilitate the acquisi- 
tion of a spectrum with sufficient signal-to noise. Through 
magnification bias, these flux limits lead to a bias on the 
probability with which GRB afterglows enter the sample as 
a function of magnification, and hence to a bias of impact 
parameter relative to the nearest galaxy. The result is that 
magnification bias concentrates observed so urce in a fiux 
limited sample around foreground objects ( [Webster et al.| 
19881, which in the context of Mg-II absorbers associated 



with galaxies should lead to an excess of absorbers per unit 
redshift in a gravitationally lensed sample. 

To quantify the potential magnitude of the multi-band 
magnification bias we specify two cumulative luminosity 
functions 'i'-,{L^) and ^A(i'A) describing the number den- 
sities of sources above fixed apparent gamma-ray Lj and 
optical La luminosities. With these the bias for observing a 



^ Gravitational lensing also lowers tlie density by magnifying the 
angular extent of the image plane relative to the source plane. 
This effect, which is usually referred to as depletion, is not dom- 
inant for steep luminosity functions. 



lensed source relative to the unlensed population is ( Wyithe 
et al.|2003| l 

1 rMmax 

B= / d^i — *T,(L^/^)*A(iA/At) 



*t,(L^)*a(La) 



dji 



(6) 

where dP/dfi is the probability distribution for magnifica- 
tion within the range ^min < fJ. < /^max- We then calculate 
the bias for the singly imaged Mg-II absorbers using the 
distribution, 



dPt 



Mg 



dfi 



for < /i < 2, (7) 



where = 1 + SerDci / R, as well as the bias for each of the 
two multiple images using the distributions 

dfj. 



(M-1)3 



for 2 < fj. < (X 



(8) 



and 



dP, 



m,2 



dfl {n + 1)3 



for < fi < oo. 



(9) 



Note that we treat the two images of a multiply imaged 
GRB separately since they have a time delay which is typ- 
ically longer than the duration of the afterglow decay time. 
To conserve flux we approximate the magnification of re- 
gions with no Mg-II absorber using a constant magnifica- 
tion (/inoMg < 1) which normalises the mean of the overall 
magnification distribution to unity. 

By combining the above expressions for the optical 
depth and magnification bias we next calculate the corre- 
sponding biased optical depths for spectra that do not con- 
tain a Mg-II absorber 



^noMg — -SnoMg 'TnoMg 



(10) 



as well as the biased optical depths for spectra that contain a 
Mg-II absorber within a spectrum of a singly imaged source, 
and within the spectra of a multiply imaged source 



Bs= da 
Jo Jo 

Bm,l ~ Bm,lT"Mg 
Sm,2 ~ -Bm,2T"Mg 



dz 



d TMg 

da dz 



BMe{z,a)TMg{z 



(11) 



In the above Bs involves integrating the differential optical 
depth over redshift and velocity dispersion because the bias 
is dependent on these quantities, whereas for multiple imag- 
ing the magnification probability distribution (and hence 
bias) are independent of these quantities. From these biased 
optical depths we are able to obtain the corresponding frac- 
tions of observed GRB afterglow spectra that do not contain 
an Mg-II absorber 

SnoMg 



-FiioMg 



(12) 



^3noMg + Bb + Bm.l + Bni,2 

as well as the fractions that contain a Mg-II absorber within 
a spectrum of a singly or multiply imaged GRB afterglow, 

Bb 



Fb = 



^noMg + Bb + Bn,,l + B 

m,2 



(13) 
(14) 



BnoMg + Bs + Bm.l + Bm,2 

The total fraction of GRB lines of sight which contain a 
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Frac. of GRB afterglow lines of sight with Mg-ll Obs. mult, image frac. GRBs with afterglows Relative numbers of Mg-ll absorbers 




Figure 1. Panels showing contours over the parameter space of a\ and imax for i) the fraction of GRB lines-of-sight with Mg-II 
absorbers (left), ii) the fraction of multiply imaged GRBs among the afterglow sample (centre) and iii) the ratio of absorbers per unit 
redshift predicted for GRB afterglow and quasar spectra (right). The model assumed a value of Rq = 165kpc, for which the fraction of 
lines-of-sight which would contain a Mg-II absorber in the absence of lensing is 0.25. In the left and right panels, the heavy contours 
correspond to the mean and 1-sigma uncertainties on the quantities -FMg and ^GRB.q respectively. 

Mgll absorber is Fms = + -Fm- In the following section we 
use the above model to calculate the fraction of lines-of-sight 
with Mg-II absorbers among the GRB afterglow sample. 



of the individual single-band luminosity functions 
f La \ 

*^,A(i7,iA) oc {L^)~°"' where La > -/^broak 

\ ^break / 

'^^,a{L^,La) oc (L^)""^ where La < ibrcak 

(15) 



3 STATISTICS OF ABSORBERS IN GRB 
AFTERGLOWS 

To calculate the fraction of GRB lines-of-sight with Mg-II 
absorbers we must specify luminosity function models for 
both the GRBs and their optical afterglows. The gamma- 
ray and optical luminosities are thought to be independent 
i Nardini et al. [2006 , 2008) , and so we specify the luminosity 
functions in these bands separately. For the GRBs we as- 
sume a shallow power-law of the form ^'^(L^) oc (L-y)""^, 
and set = 0.7 to correspond to the measured differen- 
tial slope of ~ 1.7 ( jSchaefer et al.| |2001). The luminosity 
function of optical afterglows is much less well-constrained. 
The observed distribution appears bimodal and relatively 



narrow, especially for the high luminosity branch (Nardini 
et al.|2008| |. Furthermore, many (~ 40% of long GRBs) are 
not observed to have optical counterparts. The latter fact 
suggests that the slope of the optical afterglow luminosity 
function could be steep. Therefore, for the optical afterglows 
we also assume a power-law 'I'a(^a) oc (La)~"*, but tak- 
ing into account the possibility that the slope is very steep. 
We note that for cumulative slopes steeper than qa > 2 the 
magnification bias becomes formally infinite. This would re- 
sult in an unobserved plethora of multiply-imaged GRBs. 
We therefore introduce a break in the luminosity function 
below which the slope is assumed to be very fiat so that 
intrinsically faint objects do not contribute to the magnifi- 
cation bias. Physically, this is also necessary to ensure that 
the number density of afterglows does not diverge. If we as- 
sume there are no sources fainter than Lbrcak, we can char- 
acterize the magnification required to view sources with the 

break luminosity Lbrcak as /imax = ^max = ilim,A/ibrcak- 

Thus given an observed luminosity limit I/um.A the intrinsic 
multi-band luminosity function is assumed to be the product 



3.1 Mg-II absorbers in quasar lines-of-sight 

In order to compare the statistics of GRB and quasar Mg- 
II absorbers we must also compute the fraction of quasar 
lines-of-sight that contain a Mg-II absorber (-FMg.q). This 
modelling is based on the model of GRB afterglow statistics 
described above, with the following modifications. Firstly, 
since quasars are selected in correlated (optical) bands there 
is only a single factor of magnification bias. Moreover we as- 
sume that the fiux limit in these samples corresponds to 
quasars that are at the faint end of the luminosity func- 
tion, as is appropriate for large modern quasar surveys such 
as the SDSS (we assume a cumulative slope of = 0.4; 



Groom et al. 20091. For such a shallow slope, the magnifi- 



cation bias is negligible; single-band magnification bias only 
becomes important for a cumulative slope of Oq ^2. In ad- 
dition, since quasars are not transient, the two images of a 
strongly lensed source are considered together rather than 
as separate sources. We note as a consistency check that our 
model predicts a multiply imaged lensed fraction of 0.001 
for quasars at z ~ 2 which is consistent with previous work 
(e.g. [Pindor et al.|2003| ). 

3.2 Results 

To produce a realisation of our model we therefore need 
to specify three parameters. Firstly there is the value of 
-Ro, which in our fiducial case we set to 165kpc. This value 
reproduces the observed F^^^^ = 0.28 for the fraction of 
quasar lines-of-sight with a strong Mg-II absorber (and is 
consistent with the high resolution spectroscopy of IProchter] 
|et al.|2006j which yielded 22 absorbers along 91 quasar lines- 
of-sight). We then compute models for a range of values of 
the very uncertain parameters describing the GRB afterglow 
LF, Qa and Zmax- The results are presented in the left and 
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Figure 2. Scatter plots of points in the model parameter space {a\, (max)- In the left panel we show the relation between the fraction 
of GRB lines-of-sight with Mg-II absorbers, and the fraction of multiply imaged GRBs. In the right panel we show the relationship 
between the fraction of GRB lines-of-sight with Mg-II absorbers, and the ratio of absorbers per line-of-sight between GRBs and quasars. 
In each case we show show the relationships assuming 3 values of _Rq in order to show the sensitivity to this parameter. To show the 
correspondence of the model to observed absorber statistics we have plotted in each panel contours at 60% and 10% of the peak likelihood 
given the constraints = 0.6 ± 0.15 and i^Mg% = 0-28 ± 0.05. 



central panels of Figure [T] where contours of the fraction of 
GRB lines-of-sight with Mg-II absorbers, and the fraction 
of multiply imaged GRBs among the afterglow sample are 
shown respectively. As expected there is clear degeneracy 
between the parameters /max and oa since large values of 
each parameter lead to large magnification biases and hence 
to large fractions of lines-of-sight with Mg-II absorbers as 
well as large multiple image fractions. 



To more clearly understand the relationship between 
the fraction of GRB lines-of-sight with Mg-II absorbers, 
and the fraction of multiply imaged GRBs among the af- 
terglow sample we plot their corresponding values for all 
points in the model parameter space (aA,/max) in the left 
panel of Figure [2] Figure [2] demonstrates the small amount 
of uncertainty introduced into this relationship through the 
unknown parameters (aA,Zmax), and indicates that larger 
fractions of lines-of-sight with Mg-II absorption correspond 
to larger multiple image fractions. This strong correlation 
is unsurprising, since Fms = Fa + Fm- However, note that 
both Fs and Fm depend on the lens model. In particu- 
lar, Fs falls as lensing becomes more important since if 
either Fug — ^ 1, Jm — >■ 1, or Fs — >■ the unlensed or 
singly imaged lines of sight become negligible even though 
Fmg = Fs = ^Mg.q ~ 0.3 in the absence of lensing. In Fig- 
ure [2] we also show corresponding sets of points for 2 addi- 
tional values of Ro lying within the range of uncertainty 
governed by the observed quantity -FMg!q ~ 0.28 ± 0.05 
(where the uncertainty is estimated from the number of 
strong absorbers in the high resolution spectroscopic sam- 
ple of Prochter et al.|[2006 \ in order to show the sensitivity 
to this parameter. Given the allowed values of _Ro, we can 
read off the predicted multiple image fraction from the upper 
panel of Figure [2] which shows that having Mg-II absorbers 
in ~ (60 ± 15)% of GRB afterglow lines-of-sight indicates a 
multiply lensed fraction of Fm ~ 20% — 60%. 



4 COMPARISON OF ABSORBER STATISTICS 
ALONG QUASAR AND GRB AFTERGLOW 
LINES-OF-SIGHT 

As alluded to in the introduction, a hitherto unresolved puz- 
zle has been the fact that quasar lines-of-sight show a factor 
of -RcRB.q ~ 2.1±0.6 fewer Mg-II absorbers per unit redshift 
than do GRB afterglow lines-of-sight. Gravitational lensing 



has been proposed as a possible solution to this puzzle ( For- 



ciani et al.|2007 Tejos et al. |2009J , though no quantitative 
calculation has been provided. In this section we use our 
gravitational lensing model to determine whether gravita- 
tional lensing can explain the difference in absorber density, 
and if so, for what regions of parameter space describing the 
unknown GRB afterglow luminosity function. 

In the right hand panel of Figure [T] we show contours 
of the ratio of absorbers per unit redshift predicted for 
GRB afterglow and quasar spectra in our model (-RcRB.q ~ 
FMg/FMe,q) over the parameter space of qa and /max. We 
find that ratios of -RcRB.q ~ 2 can be obtained for cumula- 
tive slopes of OA ~ 3.5 or steeper provided the power-law 
extends sufficiently far below the flux limit. As before there 
is degeneracy and to more clearly understand the relation- 
ship between the fraction of GRB lines-of-sight with Mg-II 
absorbers and the ratio of absorbers per line-of-sight be- 
tween GRBs and quasars we plot their corresponding values 
for all points in the model parameter space (aA,/max) in 
the right panel of figure [2] As before we have plotted this 
relation for three values of Rq. 



5 ABSORBER REDSHIFT DISTRIBUTION 

In addition to the statistics of the number of absorbers rel- 
ative to quasars, it is important that the model also suc- 
cessfully reproduce the observed absorber redshift distribu- 
tion. The predicted differential and cumulative redshift dis- 
tributions for a fiducial model (_Ro — 165kpc) of oa ~ 4 



and /max = 50, giving 7?grb,< 



2.4 and 



0.3, are 



shown in the left and central panels of Figure [S] respectively. 
Three model distributions are shown, with the absorber red- 
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Mg absorber redshift distributions Mg absorber redshift distributions Observed magnification distribution 
-i 1 1 1 1 1 1 1 1 I 1 1 1 1 1 1 1 1 I ' ' " I ' ' — I 1 ' — '— 




100 



Figure 3. The left and central panels show differential and cumulative absorber redshift distributions respectively. We have picked a 
fiducial model of ax = 4 and imax = 50 for this example. Three model distribution are shown corresponding to our absorber model (solid 
lines), to the distribution of gravitational lens redshifts (long dashed lines), and a distribution that assumes no lensing. The cumulative 
distributions (central panel) are compared to the observed absorber redshift distribution (thick stepped line) and the probabilities 
associated with a KS test between the models and data are listed. The observed (solid line) and intrinsic (dashed line) magnification 
distributions for the fiducial model is shown in the right hand panel. 



shift distribution strongly constrained a-priori by the path- 
lengths probed by the spectra in each case. The solid curves 
represent the absorber redshift distribution predicted by our 
model. 

In addition to the Mg-II absorber model, we show for 
comparison the distribution of gravitational (multiple im- 
age) lens redshifts (long dashed lines). This distribution is 
more peaked than predicted by the absorber model which 
includes both lensed and unlensed sources. We also show 
a model that assumes no lensing, which leads to a distri- 
bution that is more spread out in redshift. The cumulative 
distributions (central panel) are compared to the observed 
absorber redshift distribution (thick stepped line). We have 
performed a KS test between each of the three model distri- 
butions and the observed distribution. The probabilities of 
finding a distribution of absorber redshifts drawn from the 
model that are less consistent with the parent distribution 
than the observations are listed. While the data is unable 
to strongly distinguish between the different models, the hy- 
pothesis of all Mg-II absorbers lying along strongly lensed 
lines of sight is somewhat less consistent (Pks = 0.04) than 
the mixed absorber or no-lens hypothesis (Pjjs = 0.07,0.09 
respectively) . 



6 OBSERVED MAGNIFICATION 
DISTRIBUTION 

In addition to modifying the fraction of sources in a flux lim- 
ited sample that are seen as multiple images, gravitational 
lensing magnification bias also modifies the distribution of 
source magnifications relative to the intrinsic magnification 
distribution. The predicted distribution of observed magni- 
fication for the fiducial model [ax ~ 4, /max ~ 50) is shown 
in the right hand panel of Figure [3] There are 3 features of 
note in this distribution. The two peaks near magnifications 
of unity result from our assumption that a magnification 
within i? of a galaxy is given by the SIS, while beyond R 
the magnifications all have a constant value (leading to the 
delta function seen in the figure). A more sophisticated anal- 
ysis would draw on a global distribution for lensing by large 
scale structure (e.g. Fluke et al.|2002 1, however we note that 



since magnifications near unity lead to negligible magnifica- 
tion bias, this does not affect the results of our paper. 

In addition to these low magnification peaks, our model 
also predicts a prominent additional peak at high mag- 
nification. This peak differs from the generic shape of 
dP/dfi oc fi~^ for strong lenses owing to magnification bias 
which greatly increases the number of sources present in 
the flux limited sample from regions of high magnification. 
As aA,/max iucrease, the number of such highly magnified 
sources increase as well. The sharp cutoff in the distribu- 
tion beyond /i = 50 corresponds to our assumed value for 
/max. It is interesting to note that the GRB afterglow lu- 
minosity function is thought to be bi-modal ( jNardini et al.| 
[2008). Naively, the bi- modality of the magnification distri- 
bution suggests a possible lensing origin for bimodality in 
the observed GRB afterglow luminosity function. However 
when the intrinsic magnification distribution is convolved 
with the luminosity function it is easy to show that lens- 
ing cannot produce a bimodal observed luminosity function 
where the intrinsic luminosity function is monotonically de- 
creasing with luminosity (instead, a large magnification bias 
leads to a slope of oa = —2 above the flux limit). We note 
that the observed bi-modal distribution is not inconsistent 
with the lensing hypothesis, because the steep power-law 
need only be present below the flux limit in order to pro- 
duce a large magnification bias. 



7 THE NUMBER OF STRONGLY LENSED 
GRB AFTERGLOWS 

One striking prediction from our model is that if gravita- 
tional lensing is responsible for the observed difference in the 
number of Mg-II absorbers per unit redshift in GRB after- 
glow and quasar spectra, then of order 30% (~ 20% — 60%) 
of the GRBs with afterglow follow up should have been mul- 
tiply imaged. This figure is consistent with the estimate of 
60% fromlTejos et al.lpOOgf for the number of GRBs that 



must have been included in the sample owing to magnifica- 
tion in order to have produced the observed PcRB.q = 2 — 3, 
and implies that ~ 9 (~ 5 — 15) of the current sample of 26 
spectroscopic afterglows were multiply imaged. On the other 
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hand no multiply imaged GRBs have been discovered, ei- 
ther within the overall GRB sample, or within the afterglow 
subsample. In this section we therefore discuss whether cur- 
rent observational campaigns would have missed this high 
multiple image fraction among the spectroscopic afterglow 
subsample. 

Unlike multiply imaged quasars, multiple images of 
GRBs or their afterglows would be seen as repeating bursts 
since the timescale of variability is shorter than the typ- 
ical time delay between multiple images. As a result, for 
serendipitous identification of a multiply imaged GRB, the 
images would need to be discovered separately. In particu- 
lar, assuming that the location of a GRB is not monitored 
beyond the time when the afterglow is no longer visible, the 
identification of a repeating image of a lensed source would 
be via a GRB trigger. The relevant statistical question is the 
following: in the case where a GRB is strongly lensed (i.e. 
multiply imaged) with one image having been located, in 
what fraction of cases would a second image also be detected 
at a later time? Now a follow up image can only be detected 
in cases where it is the leading image that is discovered first. 
The fraction of detected images of lensed sources that are 
leading images is /icad = Sm,i/(Bm,i -I- 8^,2)- For our fidu- 
cial model we find /icad = 0.56, which is slightly greater 
than 1/2 owing to the increased effects of magnification bias 
for the leading image. Assuming monitoring coverage that 
is random with respect to the unknown time delay and that 
a GRB monitor surveys Hgrb sr at any one time, there is 
a f2GRB/47r chance that the telescope would be monitoring 
the region of sky where the lead image of the GRB was 
identified at this later time. The probability that a partic- 
ular multiply imaged GRB would be identified is therefore 
-Pm = /ieadf2GRB/47r ~ 0.09. Glvcu a Sample of Ngkb = 26 
afterglow spectra, with a Fm = 30% lens fraction, the prob- 
ability of identifying n lenses among an intrinsic sample of 
Niena = FmNQ-nB IS glvcu by the binomial theorem 



P{n) 



n[{Ni,n.-ny. 



JVi„„,-n 



(16) 



For the Swift telescope, which monitors JIgrb ~ 2 sr of the 
sky the probability of finding no lenses is therefore P(0) — 
47%, while the probabilities of having found 1, 2 or 3 lenses 
are P(l) = 37%, P(2) = 13% and P(3) = 3%. Thus we 
find it unsurprising that lensed GRBs among this sample 
have not been identified as repeating bursts, even prior to 
considerations regarding the difficulty of their identification. 



For instance, a search by Veres et al. ( 2009 \ for lensed bursts 



in the Fermi GRB data found 3 candidates with plausible 
timescales, but in the end no definitive lensing detections. 

The above considerations apply to random surveys of 
the GRB sky. On the other hand, the very large predicted 
lens fraction implies that targeted optical monitoring of the 
GRBs with multiple Mgll absorption in their afterglows in 
the months and years following the decay of the initial after- 
glow light curve would be a very efficient way of identifying 
these strongly lensed GRBs. This would provide a crucial 
test of our model. Lack of identification of lensed afterglows 
among this sample would indicate that gravitational lens- 
ing is unable to explain the difference in absorber statistics 
between quasars and afterglow lines-of-sight, and therefore 
that the afterglow LF is not steep (qa > 3) below the spec- 
troscopic flux limit. 



Finally it is worth noting that in contrast to strong Mgll 
systems with Wr 1 A, weak Mgll systems with 0.07 — 
0.3 A ^ Wi ^ 1 A have comparable incidence rates towards 
GRB and QSO fines of sight (Tejos et al.[2009l[Verg"ani et al. 
[20Q9| . The incidence of CIV also appears to be similar along 
GRB and QSO lines of sight ( [Tejos et al.|2007| ). This agrees 
well with our model under the reasonable hypothesis that 
CIV and weaker Mgll absorbers are associated with lower 
density gas which lie further away from galaxies, and as such 
are not subject to lensing bias. 



8 SUMMARY 

In recent years the rapid optical imaging and spectroscopic 
foUowup of GRBs and their afterglows has began to offer a 
new probe of the intergalactic medium out to high redshift 
which complements the more abundant quasar lines-of-sight. 
Studies of the incidence of Mg-II absorption systems along 
these lines-of-sight have revealed that although both quasars 
and GRB afterglows should provide a-priori random sight- 
lines through the intervening IGM, strong Mg-II absorbers 
are several times as likely to be found along sight-lines to 



GRBs ( |Prochter et al.|2006| [Tejos et al.|2009l|Vergani et al 
20091. Several proposals to reconcile this discrepancy have 
been put forward (see e.g. Porciani et al. 20071, but none 



have been quantitatively successful. In this paper we have 
described a simple model that associates Mg-II absorbers 
with the foreground galaxy population, and includes a 2- 
band luminosity function to describe the number densities 
of GRBs and their afterglows. We have used this model to 
estimate the effect of gravitational lensing by galaxies and 
their surrounding mass distributions on the statistics of Mg- 
II absorption. 

Our model leads to two main findings. Firstly, we show 
that the multi-band magnification bias could be very strong 
in the GRB afterglow population. Gravitational lensing can 
explain the discrepancy in the incidence of absorbers be- 
tween quasar and GRB lines-of-sight for GRB afterglow lu- 
minosity functions with cumulative slopes oa ^ 3.5. Sec- 
ondly our model makes the prediction that approximately 
20%-60% (i.e. between ~ 5 and 15) of the existing afterglow 
sample would have been multiply imaged, and hence repeat- 
ing sources. We show that despite this large lensing fraction 
it is likely that none would yet have been identified by chance 
owing to the finite sky coverage of GRB searches. As a result 
we predict that continued optical monitoring of the bright 
GRB afterglow locations in the months and years follow- 
ing the initial decay would lead to identification of lensed 
GRB afterglows. A confirmation of the lensing hypothesis 
would allow us to constrain the GRB luminosity function 
down to otherwise inaccessibly faint levels, with potential 
consequences for GRB models. 
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